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Abstract 

With  the  advent  of  nanosat  and  picosat-size 
satellites,  the  design  of  miniature,  efficient,  high  thrust- 
to-weight  (T/W)  ratio  thrust  vector  control  (TVC) 
thrusters  with  liquid  bipropellants  has  been  evolving. 
The  extremely  high  T/W  ratios  achievable  with  micro 
miniaturization  of  liquid  bipropellant  engines  are  a 
result  of  the  “cube-square-law”.  As  the  engine  is  scaled 
down  linearly,  the  propellant  flow,  and  thus  the  power, 
decreases  with  chamber  cross-sectional  area  or  the 
square  of  the  linear  size,  while  the  weight  decreases 
with  the  volume  of  the  engine  or  the  cube  of  the  linear 
size.  As  a  result,  the  power-to-weight  ratio  increases 
linearly  as  the  engine  size  is  reduced.  Only  recently,  by 
going  to  micron-scale  microelectromechanical  system 
(MEMS)  silicon  wafer  fabrication  technology,  has  the 
increase  in  T/W  ratio  at  reduced  bipropellant  engine 
scales  been  realized  over  existing  T/W  ratios  for  TVC 
thrusters.  MEMS  and  possibly  nanotechnology 
fabrication  techniques  are  the  enabling  technologies  for 
achieving  these  higher  propulsive  system  T/W  ratios.  A 
fundamental  problem  with  achieving  high  specific 
impulse  or  thrust  performance  at  these  reduced  engine 
scales  is  the  combustion  residence  time.  By  increasing 
combustion  chamber  pressure  and  adjusting  oxidizer  to 
fuel  mixture  ratio  to  fuel  rich  at  these  reduced  engine 
scales,  efficient  bipropellant  combustion  is  shown  to  be 
analytically  and  experimentally  achievable  at  engine 
characteristic  lengths  one  to  two  orders  of  magnitude 
smaller  then  characteristic  lengths  of  typical  spacecraft 
or  satellite  TVC  thrusters.  The  advantages  of  using 
reduced  scale  TVC  thrusters  in  improving  reentry 
vehicle  (RV)  targeting  accuracy  and  mission  impulse 
requirements  will  be  discussed  with  RV  and  nano  RV 
reaction  control  systems. 

I.  Introduction 

The  extremely  high  T/W  ratios  achievable 
with  micro-miniaturization  of  liquid  bipropellant 
engines  are  a  result  of  the  “cube-square-law.”1  As  the 
engine  is  scaled  down  linearly,  the  propellant  flow,  and 
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thus  the  power,  decreases  with  chamber  cross-sectional 
area  (the  square  of  the  linear  size)  while  the  weight 
decreases  with  volume  of  the  engine  (the  cube  of  the 
linear  size)  so  that  the  power-to-weight  or  T/W  ratio 
increases  linearly  as  the  engine  size  is  reduced.  Two 
different  fabrication  technologies  are  investigated  in 
pursuing  this  high  T/W  goal.  The  first  is  MEMS 
silicon  wafer  technology  being  applied  to  micro 
propulsion  research  by  the  Massachusetts  Institute  of 
Technology  (MIT).  The  second  approach2,  developed 
by  the  Air  Force  Research  Laboratory  (AFRL), 
involves  micro-machining  conical  nozzle  geometries 
with  circular  throats  down  to  approximately  50pm. 
The  AFRL  fabrication  approach  is  applicable  to  high- 
temperature  rare-earth  materials  typically  used  in 
radiation-cooled  TVC  nozzles  such  as  rhenium,  or 
tungsten.  The  primary  advantage  of  the  AFRL  micro¬ 
nozzle  fabrication  approach  is  the  capability  of 
surviving  the  corrosive  environments  of  storable 
hypergolic  propellants,  such  as  hydrazine,  used  in 
satellite  TVC  or  deep-space  propulsion  applications.  A 
proof  of  principle  micro-nozzle  of  nickel  fabrication 
with  the  AFRL  approach  is  shown  in  Figure  1 . 


Figure- 1.  AFRL  Micronozzle  Schematic.  All 
Dimensions  in  pm.  (Ketsdever,  1999) 

The  MIT  MEMS  fabrication  approach  results  in  a  high 
aspect  ratio  rectangular  bipropellant  engine  assembled 
from  silicon  wafers.  MIT  has  developed  MEMS 
fabrication  procedures  for  turbo-pumps,  regenerative- 
cooled  rocket  combustion  chambers  and  nozzles, 
propellant  injectors,  and  has  also  experimentally 
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demonstrated  these  components3'4.  The  general 
arrangement  of  the  MIT  engine  is  shown  in  Figure  2. 


Valves 


Figure-2.  MIT  MEMS  Bipropellant  Engine  General 
Arrangement  (London, 2000) 

By  modeling  the  combustion  chamber  and 
nozzle  as  a  pressure  vessel5,  the  process  of  achieving 
high  T/W  ratio  can  be  understood.  The  following 
pressure  vessel  equation  describes  how  the  thrust 
chamber  T/W  ratio  increases  with  reduction  in  scale. 

F  Pc^T^f  Cj  1 

-  OC  - —  oc - —  = - oc  -  (1) 

m  PsSAL  1}  L  L 

Once  the  chamber  pressure  reaches  a  certain  level,  and 
if  the  exit  nozzle  has  a  large  exit  to  throat  area  ratio, 
then  the  thrust  coefficient  really  is  not  changing 
dramatically,  so  the  principle  mechanism  for  chambers 
T/W  increase  is  reduction  in  chamber  characteristic 
length.  Historical  engine  weight  trends  do  not  support 
this  miniaturization  hypothesis  when  engine  chamber 
T/W  ratio  is  plotted  versus  engine  chamber  weight  as 
shown  in  Figure  3. 


Figure-3.  Bipropellant  Thrust  Chamber  T/W  Ratio 
versus  Chamber  Weight 


Figure  3  indicates  that  for  smaller  and  smaller  thrust 
chamber  weights,  the  lower  the  chamber  T/W  ratio. 
Except  for  the  MIT  MEMS  engine  and  the  Marquardt 
LEAP  thruster,  this  trend  seems  counter-intuitive  to  the 
analytical  argument  put  forth  in  Equation  1.  The 
reason  for  this  trend  reversal  can  be  explained  by 
plotting  chamber  T/W  ratio  versus  the  Bayt 


characteristic  length  parameter5, 


as  shown  in 


Figure  4. 
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Figure-4.  Bipropellant  Thrust  Chamber  T/W  Ratio 
versus  Cf/L* 


Figure  4  indicates  that  very  few  rocket  engines  have 
really  small  characteristic  lengths  representative  of  the 

MIT  MEMS  rocket  engine  ( L  ~  .25  feet )  .  Most 


bipropellant  engines  have  a 


of  approximately 


0.55  with  a  characteristic  length  around  3  to  4  feet.  In 
order  to  achieve  truly  high  T/W  ratios,  bipropellant 
engines  of  very  short  characteristic  length  have  to  be 
developed  to  take  advantage  of  the  pressure  vessel 
scaling  equation  (Eq.  1).  Either  the  MIT  MEMS  or  the 
AFRL  metal  micro-nozzle  fabrication  approach 
achieve  the  high  T/W  ratios  predicted  by  Equation  1. 
MIT  has  experimentally  proven  that  bipropellant  rocket 
T/W  ratio  increases  at  reduced  engine  characteristic 
length  (Figure  3)  but  not  as  high  as  originally 
proposed3  (T/W  equal  to  1200).  MIT  would  have 
experimentally  demonstrated  a  much  higher  T/W  ratio 
by  going  to  thinner  gauge  silicon  wafer  material  then 
used  on  their  prototype  engine  which  was  designed  to 
operate  at  .125  atm.  The  experimental  tests,  only 
achieved  a  combustion  chamber  pressure  of  12  atm. 


To  take  advantage  of  the  enhanced  T/W  ratio 
at  reduced  engine  scales,  several  issues  that  diminish 
performance  at  reduced  engine  scales  have  to  be 
addressed.  The  first  issue  is  the  reduced  chamber 
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residence  time  at  reduced  engine  scales  requiring 
shorter  and  shorter  liquid  bipropellant  combustion 
residence  times.  The  chamber  residence  time  is  given 
by  the  following  equation6 


(2) 


AV  =g0IspIn^L-frgcosM 


m 


find 


C d^-ref  T  P atmV2 

-  - - - dt  -  — (1  -  cos  ccc  )dt 


6 rF ■ 


m 


m 


(4  b) 


Reference  6  further  notes  that 


Pcj 


m 


is  independent 


of  chamber  pressure  such  that  chamber  residence  time 
is  independent  of  chamber  pressure,  but  combustion 
residence  time  is  not.  By  increasing  the  chamber 
combustion  pressure,  and  adjusting  the  oxidizer  to  fuel 
ratio  to  fuel  rich  as  described  in  the  theoretical  section, 
the  combustion  residence  time  can  be  shortened  to  the 
desired  values  below  the  chamber  residence  time. 


Another  important  factor  at  reduced  engine 
length  scales,  is  the  high  heat  loss  in  relation  to  heat 
generated  during  combustion7.  The  heat  loss  in 
relation  to  heat  generated  is  inversely  proportional  to 
the  engine  characteristic  length  according  to  the 
following  relationship. 


0) 


where  E *  is  the  heat  loss  in  relation  to  the  heat 

generated,  E .  At  the  high  power  to  weight  ratios 
associated  with  the  small  engine  characteristic  lengths 
utilizing  MEMS  fabrication  techniques,  regenerative 
cooling  would  be  required  to  recover  the  high  heat 
losses  during  combustion1.  Therefore,  with 
regenerative  cooling,  high  specific  impulse  operation 
can  be  achieved  with  very  small  engine  characteristic 
lengths,  lengths  on  the  order  of  0.1  to  1  inch. 


The  next  performance  reducing  phenomena  at 
reduced  length  scales  is  the  aerodynamic  losses 
associated  with  small  scale  rocket  vehicles8.  To 
understand  the  effect  of  scale  upon  the  effective 
specific  impulse,  the  rocket  delta-V  equation  is  used 
with  losses8  ,or 


AK  delivered  gravity  drag  ^  ^ steering 

(4.  a) 


For  very  small  rocket  vehicles,  the  vehicle  mass,  m,  is 
small  in  relation  to  the  cross-sectional  drag 
area ,CdArej ,  This  reduction  in  vehicle  ballistic 


coefficient, 


m/ 


' CdAref 


causes  AV  losses  in  vehicle 


performance  for  very  small  vehicles.  Figure  5  shows 
the  effect  of  air  drag  upon  launch  vehicle  performance. 
As  the  launch  vehicle  gross  takeoff  weight  is  reduced, 
the  drag  loss  portion  of  total  AV  increases,  and  the 
payload  mass  fraction  is  reduced.  Figure  5  indicates 
that  the  larger  the  launch  vehicle  ballistic 

coefficient,  m/r  A  ,  the  larger  the  vehicle  payload 
/  C DAref 

mass  fraction.  This  implies  that  in  order  to  take 
advantage  of  the  high  T/W  ratio  performance 
enhancement  of  small  scale  rocket  engines,  they  should 
be  clustered  into  larger  arrays  propelling  larger 
vehicles1.  In  this  way,  the  reduction  in  propulsive 
system  weight  in  relation  to  the  overall  vehicle  weight, 
would  be  realized  without  a  corresponding  loss  in 
payload  mass  fraction  due  to  the  enhanced  effects  of  air 
drag. 


-  Payload  Fraction  of  Vehicle  Mass  -  ■-  Drag  Loss  Fraction  of  Total  Velocity  Increment 
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Figure-5.  Effects  of  Air  Drag  on  Launch  Vehicle 
Performance  (770  km  Orbit/90  deg.  Inclination/Ground 
Launched)  (Francis,  1999) 

A  beneficial  implication  of  the  delta-V 
equation  for  small  vehicles  is  during  atmospheric 
reentry.  While  the  delta-V  equation  (Eq.  4.b)  states 
that  payload  mass  fraction  of  a  reduced-scale  launch 
vehicle  is  significantly  reduced,  this  is  not  the  case  for 
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a  reduced-scale  weapon  payload  such  as  a  nano  RV.  If 
you  eliminate  the  thrust  terms  of  the  delta-V  equation 
(Eq.  4b),  and  translate  the  coordinate  system  from  a 
vertical  launch  system  to  a  horizontal  reentry  vehicle 
coordinate  system,  the  equation  reduces  to  the 
deceleration  equation  for  a  RV,  or 


-dV 

dt 


-g  sin  6  + 


pv2 

m  2 


(4.c) 


For  a  reduced-scale  RV  ,  such  as  a  nano  RV  ,  the 
vehicle  decelerates  more  quickly  during  atmospheric 
reentry  due  to  the  lower  ballistic  coefficient,  mfCDArej . 

This  results  in  the  nano  RV  decelerating  below 
plasma  shielding  speed  more  quickly  allowing  the 
GPS/INS  guidance  system  to  acquire  targeting 
coordinates  with  more  time  before  impact.  As  a  result 
micro-miniaturization  helps  you  achieve  targeting 
objectives,  either  thru  microthruster  TVC  enhancing 
pointing  accuracy,  or  through  the  enhanced 
deceleration  of  reduced-sized  RVs.  The  advantages  of 
using  reduced-scale  TVC  thrusters  in  improving  RV 
targeting  accuracy  and  mission  impulse  requirements 
will  be  discussed  with  RV  and  nano  RV  reaction 
control  systems. 

IL  Theoretical  Analysis 


11.A.  Microthruster  Combustion  Residence  Time 
Analysis 


In  the  interest  of  finding  how  small  a 
bipropellant  engine  can  be  made  and  still  achieve 
respectable  specific  impulse  performance,  a  more 
detailed  analysis  of  the  combustion  residence  time  as  a 
function  of  chamber  pressure  and  oxidizer  to  fuel  (O/F) 
ratio  was  made.  The  combustion  residence  time 
modeling  that  AFRL  did  involved  data  matching  more 
detailed  theoretical  combustion  residence  time  analysis 
done  by  MIT9.  MIT  used  two  different  computational 
methods  described  in  References  10,  11,  and  12.  A 
more  detailed  discussion  of  the  MIT  MEMS 
bipropellant  engine  design  condition  be  found  in 
Reference  3  and  was  for  a  propellant  mixture  of 
ethanol  and  oxygen  at  a  125-atm  chamber  pressure. 
The  MIT  combustion  residence  time  analysis  results 
are  shown  in  Figure  6  versus  combustion  temperature. 


Figure-6.  Chemical  Kinetics  Residence  Time  for 
Ethanol  and  Oxygen  (O.M.  Al-Midani,  1998) 

The  interesting  result  of  the  MIT  analysis  is  that  the 
shortest  combustion  residence  time  occurs  at  the  lowest 
combustion  temperature  near  the  fuel-rich  blowout 
limit.  This  occurs  at  a  O/F  ratio  of  1.15  at  a 
combustion  temperature  of  2700  degrees  Kelvin  as 
shown  in  Figure  7. 


Figure-7.  Combustion  Chamber  Equilibrium 
T emperatures  for  Different  O/F  Ratios  for  Ethanol  and 
Oxygen  (O.M.  Al-Midani,  1 998) 

At  the  stoichiometric  O/F  mixture  of  2.2,  where  the 
highest  obtainable  combustion  temperature  occurs,  the 
combustion  time  is  orders  of  magnitude  larger  than  the 
fuel-rich  cases,  and  stoichiometric  combustion  is  not 
considered  possible  for  the  short  0.1 -msec  chamber 
residence  time  of  the  MIT  MEMS  bipropellant  engine. 
This  was  eventually  proven  by  MIT  experimentally 
when  their  MEMS  engine  only  operated  over  the 
narrow  fuel-rich  range  shown  in  Figure  10. 

II.B,  Bimolecular  Chemical  Kinetics  Combustion 
Model 
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From  the  MIT  analysis  results  in  Figure  6,  the 
combustion  residence  time  appears  to  be  proportional 
to  the  Arrhenius  activation  energy,  which  is  reduced  at 
the  lower  fuel-rich  combustion  temperatures  resulting 
in  reduced  combustion  residence  time.  This  forms  the 
basis  of  data  matching  the  MIT  combustion  analysis9 
with  a  more  general  bimolecular  chemical  kinetics 
combustion  model  utilizing  Arrhenius  activation 
energy  modeling  to  calculate  combustion  residence 
time13.  The  second-order  chemical  kinetics  rate 
equation  is  defined  as 


-dN 


dt 


i-  =  KAl-Nt-N, 


where 


~dNA 

^  _  dt 

AB  7T~ 

na-nb 


(5.a) 


(5  b) 


N A  and  N B  are  respectively  the  number  of  fuel  and 

oxidizer  molecules  per  unit  volume,  and  KAB  is  the 
second-order  combustion  rate  constant  for  unlike  fuel 
and  oxidizer  molecules.  For  combustion  of  fuel  and 
oxidizer  molecules  at  a  given  O/F  ratio,  the  following 
stoichiometric  equation  is  used. 

A  +  bB  =  eE  +  fF  (6) 

For  the  above  general  reaction  under  steady  state,  the 
rate  of  combustion  of  fuel  to  oxidizer  molecules  is 


-Ma 

dt 


l  dN  B  _dx  _ 
b  dt  dt 


■NA-NB  (7) 


X  is  defined  as  the  number  of  molecules  of  combustion 
products  created  per  unit  volume.  Integration  of  the 
above  rate  equation  yields 


f  Nb.  ~  Nb ^ 


(8) 


dx 

solving  for  — ,  we  get  the  following  equation  for 

dt 

combustion  products  generated. 


dx 

~dt 


=  KAB-b  •  (Na>  -  X )2  =  K2B  ■  {nAo  -  xf  (9) 


Equation  5. a  for  unlike  molecules,  N A  and  N B  has 
been  reduced  to  a  kinetics  rate  equation  for  like  N A 

molecules,  where  K2B  is  the  second-order  combustion 
rate  constant  for  like  fuel  and  oxidizer  molecules.  The 
solution  to  Equation  9  is14 

m  NA„-b-KAII->  n 

Ha.  (1  +  NA.-bKABt) 


where 

N^-bKAB=t-  (M) 

and  T  is  defined  as  a  combustion  time  constant.  For 
complete  combustion,  we  assume  the  ratio  of 
combustion  products  concentration  to  the  initial  fuel 
concentration  is  complete  at  85  percent  or 


or 


comb 


=  6r  =  - 


NA-b-KAB 


(12.a) 


(12  .b) 


To  calculate  the  combustion  rate 
constant,  KAB  ,  the  kinetic  theory  of  gases  is  used. 
The  number  of  bimolecular  collisions,  ZAB ,  between 

unlike  molecules  A  and  B  per  unit  volume  per  unit 
time  is  given  by  the  following  expression 

Zab  =  ^-(Tab2  -Vr- Na-Nb  (13) 


Na0  and  N B0  are  the  initial  number  of  fuel  and 

oxidizer  molecules  per  unit  volume  prior  to 
combustion.  By  rearranging  equation  5. a,  7,  and  8,  and 


where  the  average  molecular  diameter ,  (JAB ,  is 


Gab  — 


<Ja  +  <J  b 
2 


(14) 
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The  mean  relative  molecular  speed,  Vr ,  of  type  A 
and  B  molecules  is 


"8 -Kb'TcombV 

J 


(15) 


and  the  reduced  molecular  mass,  JU  ,  is 


M  = 


ma  •  mb 
(i ma  +  mb ) 


(16) 


The  average  molecular  diameter,  (jAB ,  can  be 
calculated  from  molecular  diameter  tables  in  Reference 


15. 

Since  the  number  of  molecules  which  react 
per  unit  time,  is  much  less  than  those  entering  into 
collision,  an  Arrhenius  activation  energy  function  is 
derived  to  calculate  the  number  of  reactant  molecules 
that  do  react.  This  expression  is  defined  as 


%AB 


—  ^  ,  q  ^ mix  ^ 'comb 


(17) 


where  Z'  is  the  number  of  effective  bimolecular 
collisions  per  unit  volume  per  unit  time,  and  E0  is  the 
activation  energy  threshold  that  the  colliding  molecules 
have  to  exceed  to  achieve  reaction.  R  ,  is  the 


universal  gas  constant  of  the  reactant  gas  mixture 
and  Tcomh  is  the  combustion  temperature,  S  is  the 

steric  factor,  a  collision  impact  orientation  factor  for 
the  number  of  molecules  that  have  the  proper 
orientation  for  reaction.  For  the  purposes  of  this 
theoretical  investigation  where  we  are  just  data 
matching  MIT  combustion  data,  S  is  assumed  to  be 
one.  The  combustion  activation  energy  is  chosen  to 
match  both  the  theoretical  data  shown  in  Figure  6,  and 
MIT  experimental  data  obtained  from  Reference  3  and 
6.  Finally,  by  combining  results  from  equation  5.b,  13 
and  17,  we  get  the  following  expression  for  the 
bimolecular  combustion  rate  constant  or 

£ 

AB  Na-Nb  Na-Nb  (18) 


■■X-<T2AB-Vr 


S-e1 


To  calculate  the  time  constant,  'V comb ,  in  equation  12.b, 

and  to  calculate  the  mixture  gas  constant,  Rmix ,  the 

oxidizer  stoichiometric  coefficient,  b,  in  Equation  6  is 
calculated  with  the  following  equation 


(£1 

f  M  fuel  \ 

UJ 

\  M  oxidizer  j 

(19) 


where  O/F  is  the  oxidizer  to  fuel  ratio  by  weight  and 
Mfuel  and  Moxidizer  are  the  molar  molecular  weight 

of  the  fuel  and  oxidizer.  Once  the  stoichiometric 
coefficient  has  been  calculated,  then  the  following 
expressions16  are  used  for  calculating  Rmix  .  The  mole 


fractions,  Xi ,  for  the  fuel  and  oxidizer  can  be  shown  to 


be 


**  1  +  * 


and 


where  P  is  the  combustion  pressure,  and  Pa  and  Ph 
are  the  fuel  and  oxidizer  partial  pressures.  With  the 
partial  pressures  Pa  and  Pb  from  Equation  21,  the  fuel 

and  oxidizer  freestream  densities  prior  to  combustion 
can  be  determined  with  the  perfect  gas  law,  or 


Pa 


R  •  T 

*  *  comb 


■>Pb  = 


R-B  'TCOmb 


(22) 


The  fuel  and  oxidizer  gas  constants  are  calculated  with 


R.  = -  and 

ma 


R-b 


R 

Mb 


(23) 


where  R  is  the  universal  gas  constant  and  M A  and 

mb  are  the  fuel  and  oxidizer  molar  weight. 

Finally,  the  gas  constant  of  the  combustion  gas  mixture 
is  calculated  with 


P^Ra^_Pa_Rb 
(Pa  +  Pb) 


(24) 
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Thru  an  iterative  process  with  equations  12.b  to  24,  the 
chemical  kinetics  residence  time  data  derived  in  Figure 
6  by  MIT  is  data  matched  for  oxygen  and  ethanol  over 
an  oxidizer  to  fuel  ratio  range  of  1.1  to  1.3.  By  going 
through  this  procedure,  we  obtain  the  range  of 
activation  energy  versus  combustion  temperature 
shown  in  Figure  8  for  ethanol  and  oxygen. 


Figure-8,  Activation  Energy  versus  Combustion 
Temperature  for  Ethanol/Oxygen 


V  _  tresjn 

At  Pc'  At 


(26.a) 


By  substituting 

m  =  pT  *  At  *  FY  =  Pt  *  At  *  cir 

and  Tcomb  from  equation  25,  we  finally  get 


Tcomb  *  Pf  * 
Pc 


(26.b) 


By  assuming  that  the  activation  energy  per  mole,  E0 , 

is  independent  of  chamber  pressure  at  the  chosen  O/F 
ratio  of  1.2,  then  the  combustion  chamber  residence 
time  can  be  calculated  as  a  function  of  chamber 
pressure  with  equation  12.b  or 


6 

NAo'b-KAB 


Pt'  aT 

Pc 


(27) 


1I.C.  Determining  Minimum  Chamber  Size  with  MIT 
Combustion  Analysis  and  Bimolecular  Chemical 
Kinetics  Model 

The  procedure  to  calculating  the  shortest 
combustion  residence  time  is  selecting  a  fuel -rich  O/F 
ratio  that  minimizes  the  combustion  temperature  or 
activation  energy  in  Figure  8,  and  which  is  still  above 
the  fuel-rich  blow-out  limit  calculated  by  MIT  in 
Figures  6  and  7.  For  this  example,  an  O/F  ratio  of  1.2 
was  chosen  just  above  the  blow-out  limit  of  1.15  from 
the  MIT  data  in  Figure  6  and  7.  This  gives  an 

activation  energy  ,  E0 ,  equal  to  71,500  Joules/mole 

and  will  be  the  basis  of  a  study  to  determine  the 
minimum  combustion  residence  time  as  the  chamber 
characteristic  time  is  reduced.  This  will  be  presented 
as  the  minimum  chamber  pressure  required  for 
complete  combustion  as  engine  characteristic  length  is 
reduced.  This  procedure  starts  with  chamber  residence 
time  (Eq.  2)  being  set  equal  to  or  greater  than  the 
combustion  residence  time  or 

7*  res  —  ‘Tcomb  (2^) 

From  equation  2  and  the  definition  for  characteristic 
length17,  Lc ,  the  chamber  characteristic  length  can  be 
defined  as 


From  equation  27,  the  chamber  characteristic  length, 
Lc ,  for  complete  combustion  is  inversely  proportional 
to  the  initial  concentration  of  fuel  for  a  constant  value 
0fKAB  ,  or  inversely  proportional  to  chamber  pressure. 

ILD.  Determinine  Minimum  Chamber  Size  utilizing 
the  MIT  Binropellant  Experimental  Results  and  the 
Bimolecular  Chemical  Kinetics  Model 

MIT  experimentally  demonstrated  a  MEMS 
bipropellant  engine3  similar  in  size  to  the  one  shown  in 
Figure  2.  The  experimental  results  were  for  gas  to  gas 
propellant  mixing  after  injection  into  the  MEMS 
combustion  chamber  with  methane  and  oxygen  as  the 
propellants  instead  of  ethanol  and  oxygen,  the 
propellants  in  the  analytical  combustion  study  9.  The 
procedure  to  matching  the  MIT  experimental  data 
involves  determining  the  chamber  pressure  where 
exhaust  velocity  is  maximized,  determining  the  O/F 
ratio  for  these  given  conditions,  and  then  determining 
the  Arrhenius  activation  energy  required  to  match  the 
combustion  residence  time  which  is  set  equal  to  the 
chamber  residence  time  of  the  MIT  MEMS  engine. 
The  date  matching  procedure  is  similar  to  the  data 
matching  of  the  MIT  bipropellant  chemical  kinetics 
analysis  results  shown  in  Figure  6,  but  is  done 
differently  depending  upon  the  available  experimental 
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data.  To  start  the  procedure,  the  chamber  residence 
time  (Eq.  2)  for  the  MIT  MEMS  bipropellant  engine  is 
a  fixed  quantity  independent  of  chamber  pressure  as 
discussed  in  the  introduction  and  is  equal  to  0.1  msec. 
The  combustion  residence  time  is  inversely 
proportional  to  chamber  pressure  represented  by  the 


initial  fuel  concentration,  N A  ,  as  determined  by 
equation  12.b.  By  examining  experimental  data  for 


characteristic  exhaust  velocity,  C* ,  as  a  function  of 

chamber  pressure,  and  determining  the  pressure  where 
* 

c  reaches  a  maximum,  the  combustion  residence  time 
is  then  assumed  to  be  equal  to  the  chamber  residence 
time  at  this  set  of  test  conditions.  This  assumption  is 

valid  because  c  is  a  slowly  varying  function  of 
combustion  temperature  or 


A.+] 

y+i]2{r-i) 

~2y 


(28.a) 


and  is  calculated  experimentally  with 


*  PC' At 

c  = — : — 


(28,b) 


m 

These  characteristic  velocity  results  are  shown  in 


Figure  9  for  the  MIT  microrocket  and  indicates  that  C 
reaches  a  plateau  around  7  atm. 
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Figure-9.  Experimental  Measurement  of  Characteristic 
Velocity  versus  Chamber  Pressure  for  MIT 
Microrocket  (London,2001) 

The  next  quantity  required  for  data  matching  the 
experimental  combustion  results  with  the  bimolecular 
kinetics  model  is  the  O/F  ratio  for  the  test  in  question. 
This  is  shown  in  Figure  10  and  turns  out  to  be 
approximately  2.6. 


8  i 
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Figure-10.  MIT  Microrocket  Ignition  Results  for 
Methane  and  Oxygen  (London,2001) 


Notice  that  the  MIT  microrocket  has  a  narrow  O/F 
operating  range  of  2.5  to  2.7,  well  below  the 
stoichiometric  value  of  4.0  for  oxygen  and  methane. 
This  is  just  what  you  would  expect  from  the  chemical 
kinetics  residence  time  modeling  done  by  MIT9  and 
shown  in  Figure  6.  At  the  higher  temperatures 
associated  with  stoichiometric  combustion,  the 
combustion  residence  time  is  much  larger  then  the 
chamber  residence  time.  In  fact,  during  the  MIT  tests, 
the  fuel-lean  blow-out  limit  was  reached  at  an  O/F  ratio 
of  3.0  as  shown  in  Figure  10.  In  order  to  data  match 
the  MIT  experimental  combustion  data  with  methane 
and  oxygen,  an  Arrhenius  activation  energy  of  91,500 
Joules/mole  was  required.  The  final  results  for 
minimum  chamber  characteristic  length  versus 
chamber  pressure  for  complete  combustion  from  data 
matching  both  the  MIT  analytical  chemical  kinetics 
residence  model9,  and  the  MIT  experimental 
microrocket  data3  ,  is  shown  in  Figure  13. 


Figure-1 1.  Minimum  Chamber  Characteristic  Length 
versus  Chamber  Pressure  for  Complete  Combustion 

The  resulting  specific  impulse  from  TEP  analysis12  was 
345  seconds  for  methane  and  oxygen  and  301  seconds 
for  ethanol  and  oxygen  at  the  chamber  pressures 
required  for  complete  combustion  at  these  reduced 
engine  characteristic  lengths.  The  smaller  the 
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combustion  chamber,  the  higher  the  chamber  pressure 
required  to  achieve  complete,  efficient  combustion. 
These  results  only  apply  near  the  fuel-rich  blow-out 
limit.  At  a  stoichiometric  O/F  ratio,  the  chamber 
characteristic  length  would  be  orders  of  magnitude 
larger  for  complete  combustion. 

III.  Application  of  High  Thrust  to  Weight  Ratio 
Binropellant  Thrust  Vector  Control  to  Reentry 
Vehicles  Thru  Micro-miniaturization 

We  have  chosen  the  simplified  reaction 
control  system  (RCS)  analytical  modeling  in  Reference 
18  for  evaluating  microthruster  TVC  for  two  different 
weight  class  RVs.  The  large  scale  RV,  shown  in 
Figure  12,  is  a  subscale  wind  tunnel  model  of  a  flapped 
R V  configuration  known  as  MaR  V  . 


than  typical  RVs.  Typical  pitch  moment  of  inertias  and 
lengths  were  used  for  an  RV  and  nano  RV  in 
calculating  RCS  propellant  consumption  rate. 

The  RCS  analytical  model  starts  with 
determining  the  minimum  RCS  thrust  requirements. 
These  thrust  levels  are  normally  determined  by  the 
maximum  allowable  angular  rates  during  certain  RV 
functions  such  as  star  pointing,  or  GPS  reception 
during  exospheric  coast  after  fourth-stage  bus 
deployment  prior  to  reentry  or  targeting  maneuvers  that 
require  a  given  orientation  or  attitude.  Because  of  the 
impossibility  of  imparting  an  infinitesimal  small  thrust 
amplitude  to  the  RV,  it  is  virtually  impossible  to  obtain 
a  given  attitude  with  zero  angular  rate.  As  a  result,  the 
RV  is  kept  at  a  given  attitude  within  a  dead-band 
region  established  by  a  limit-cycle  amplitude  shown  in 
Figure  14. 


TOP  VIEW 

Figure-12.  Maneuvering  RV  Details 


To  make  this  RV  controllable  with  reaction  control 
jets,  the  flaps  are  replaced  with  jet  rows.  For  RCS 
thrusters  of  very  small  size,  the  jet  rows  would  turn 
into  thruster  arrays  which  would  cover  the  conic  slice 
flat  sections  near  the  RV  frustrum.  The  second  RV  is  a 
much  smaller  RV  known  as  a  nano  RV.  A  nano  RV  is 
shown  in  Figure  1 3  during  full-scale  ablation 
aerodynamics  testing  in  the  NASA  Ames  20-mW  arc 
jet19.  Nano  RVs  are  1  to  2  orders  of  magnitude  smaller 


Figure-13.  Nano  Reentry  Vehicle  (Haldeman,  1986) 


t 

dead  band 

_ L_ 


Figure- 14.  Reentry  Vehicle  Dead  band  Motion 
(Pohl,  1965) 

To  minimize  RCS  fuel  consumption,  the  aerodynamic 
center  of  pressure  is  positioned  at  the  center  of  gravity. 
The  RV  is  not  aerodynamically  stabilized,  but  is 
stabilized  with  reaction  control  jets.  The  dead-band 
angular  tolerance  for  attitude  keeping  is  established  by 
the  attitude  accuracy  requirements. 

To  determine  the  RV  RCS  propellant 
consumption  rate  during  limit-cycle  operation,  the 
procedures  described  in  Reference  1 8  are  used.  If  the 
limits  of  the  dead  band  are  ±  (3  ,  then  the  time  history 
shown  in  Figure  14  results  for  the  RV  pitching  motion. 
The  propellant  consumption  per  second  of  operation, 
W ’  ,  for  the  limit-cycle  pitching  motion  can  be  shown 

to  be 


( Y\k 

4-Jx-Lp-P 


(29) 


where  L  is  the  reaction  jet  moment  arm,  1  is  the 
impulse  bit  (F.dt)  provided  by  the  reaction  control  jet, 

and  Jx  is  the  reentry  vehicle  polar  moment  of  inertial. 
The  attitude  keeping  propellant  consumption  rate,  W p  , 
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is  a  strong,  second-order  function  of  the  minimum 
impulse  bit,  I.  The  impulse  bit  wants  to  be  as  small  as 
possible,  especially  for  very  small  dead-hand  angular 
tolerances. 

As  an  example,  the  attitude-keeping  propellant 
consumption  for  a  reentry  vehicle  reaction  control 
system  was  calculated  during  exospheric  coast  after 
bus  deployment  and  during  reentry.  As  shown  in 
Figure  15,  by  placing  the  center  of  pressure  at  the 
center  of  gravity,  the  RCS  propellant  consumption  is 
minimized. 


Figure-15.  Reentry  Vehicle  Reaction  Control  System 
Forces  and  Moments 

This  results  in  the  aerodynamic  pitching  moment,  M p  , 
going  to  zero. 


Mp=Cna.j3- A-q-sm-l  -»  O  (30.a) 

Also,  the  aerodynamic  pitch  damping,  M  d ,  can  also 
be  neglected  as  follows 


Md=Cmq-^-q-d2  =  0  (30.b) 

The  pitch  damping  coefficient,  C  can  be  shown  to 
be  small20. 

Cmq  =  Cna.2.{smf  =  0  (30  c) 

fid 

Also,  the  induced  angle  of  attack,-^-,  during  dead¬ 
band  pitching  motion  is  very  small  at  reentry  speeds. 
Essentially,  the  reentry  vehicle  has  negligible 
aerodynamic  damping,  and  aerodynamic  pitching 
moment,  or  the  propellant  consumption  rate  for  the 
RCS  during  both  the  exospheric  and  reentry  portion  of 
the  RV  trajectory  can  be  calculated  with  equation  29. 


Another  interesting  feature  about  reaction 
control  jets  is  their  interaction  with  surrounding  air 
flow  during  reentry  as  illustrated  in  Figure  15.  This 
control  jet  interaction  increases  the  effectiveness  of  the 
control  jet21.  This  increased  effectiveness  is  usually 
referred  to  as  the  thrust  amplification  factor  and  is 
defined  as 


Fi  +  Thrust 
Thrust 


(31) 


Some  data  for  kf  as  a  function  of  temperature  and 

molecular  weight  of  exhaust  products  is  shown  in 
Figure  16. 


Figure-16.  Dependence  of  the  Jet  Thrust  Amplification 
Factor,  Kf,  on  Temperature  and  Molecular  Weight 

By  assuming  a  60-percent  increase  in  reaction  jet  thrust 
(kj  equal  to  2.6),  then  the  effective  Isp  will  be  about 

60-percent  higher  due  to  jet  interactions  with  the 
freestream  during  reentry. 

An  estimate  of  the  minimum  impulse  bit,  1, 
has  to  be  made  in  order  to  calculate  the  limit-cycle  fuel 

consumption  rate,  W  .  By  photographically  sub¬ 
scaling  the  MIT  MEMS  bipropellant  engine  (Figure  2), 
and  utilizing  the  minimum  chamber  pressure  for 
complete  combustion  versus  chamber  characteristic 
length  results  in  Figure  1 1,  an  estimate  of  the  minimum' 
obtainable  thrust  as  a  function  of  engine  characteristic 
length  can  be  made.  These  microthruster  results  are 
shown  in  Figure  17  for  methane/oxygen  and 
enthanol/oxygen  propellants. 
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Figure-17.  Minimum  Obtainable  Microjet  Thrust 
versus  Chamber  Characteristic  Length 

It  is  interesting  to  note,  that  because  the  minimum 
chamber  pressure  for  complete  combustion  versus 
chamber  characteristic  length  is  utilized  in  the 
calculation,  the  minimum  obtainable  thrust  versus 
chamber  characteristic  length  results  from  the 
calculation.  I ,  the  minimum  obtainable  impulse  bit,  is 
calculated  with  the  minimum  obtainable  thrust  as 
follows. 

I  =  |f  Fdt  (32) 

An  estimate  of  2  milliseconds  for  valve  opening  time, 
AT ,  was  made22.  Valve  leakage  rates  should  not  be  a 
major  concern  for  RVs  because  the  operational  mission 
is  usually  20  minutes  or  less  from  the  fourth-stage  bus 
dispense  to  impact.  Satellite  station  keeping  with  much 
longer  operational  missions  lasting  months  to  years, 
would  have  to  address  valve  leakage  issues. 

The  final  quantity  required  before  calculating 
the  dead-band  limit-cycle  propellant  consumption  rate, 
is  the  dead-band  angular  target  pointing  tolerance,  ft , 
in  Figure  14.  Due  to  RV  ablation  shape  change  aero 
uncertainties,  this  dead-band  target  pointing  tolerance, 
±  ,  was  set  to  ±  0.1  degree  23,24. 

Finally,  the  mission  impulse  requirements  for 
the  dead-band  limit-cycle  is  calculated  as  follows 


RV  bus.  The  trajectory  time,  AT ,  is  for  a  typical 
exospheric  suborbital  ballistic  trajectory,  and  a  typical 
atmospheric  reentry  to  impact.  The  final  mission 
impulse  energy  requirements  are  shown  in  Figure  18 
for  the  dead-band  pointing  limit-cycle  as  a  function  of 
microthrust  engine  characteristic  length. 


Figure-18.  RV  RCS  Mission  Impulse  Requirements  for 
Dead-Band  Targeting  Motion 

The  dead-band  pointing  tolerance  energy  consumption 
is  quite  small,  and  in  reality,  any  deviation  from  the 
target  flight  path  would  require  significantly  more  RCS 
total  impulse  then  shown  due  to  ablation  shape  change, 
and  changes  in  aerodynamic  center  of  pressure  during 
reentry.  What  this  mission  impulse  requirements 
analysis  in  Figure  1 8  does  show  is  for  smaller  values  of 
the  minimum  impulse  bit,  I,  the  smaller  the  mission 
impulse  requirements  to  maintain  a  given  dead-band 
pointing  tolerance.  Also,  the  nano  RV  has  higher 
mission  impulse  requirements  than  the  RV  for  a  given 
thruster  characteristic  length  because  the  microthruster 
is  initially  oversized.  Basically,  the  smaller  the  RV, 
the  smaller  the  minimum  impulse  bit  required  to 
maintain  a  given  pointing  angular  tolerance.  For 
tactical  warhead  RV  applications,  the  RV  has  to  have 
central  error  probable  (CEP)  of  35  feet  to  be 
effective24.  Having  RCS  thrusters  with  a  very  small 
minimum  impulse  bit  goes  a  long  way  to  achieving  this 
accuracy  requirements. 

IV,  Conclusions 


T  —  Y  W  AT  7  ni)  Small-scale  flight  vehicles  make  outstanding 

T  —  L  p*  /  *  Spt  v  ;  missile  payloads  such  as  a  MEMS  satellite,  or  nano  RV 

where  the  enhanced  drag  characteristics  of  a  reduced 
1SP  is  varied  for  reentry  with  jet  interaction  mass  vehicle  is  an  advantage  during  reentry.  Using 
'  MEMS  technology  to  develop  a  reduced  scale  launch 

argumentation,  or  for  vacuum  thrust  performance  vehjcle  is  counter  productive  because  of  the  reduced 
during  the  exospheric  dispense  from  the  fourth-stage  payload  mass  fraction  resulting  from  increased 
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aerodynamic  losses  at  reduced  scales.  The  enhanced 
drag  characteristic  of  nano  RVs  during  reentry  is 
beneficial,  because  they  will  decelerate  below  the 
plasma  shielding  velocity  more  quickly  allowing  more 
time  for  target  acquisition  and  seeking.  This  also 
implies  that  in  order  to  take  advantage  of  the  high  T/W 
ratio  enhancement  of  microthruster  engines,  they 
should  be  clustered  into  larger  arrays  propelling  larger 
vehicles.  In  this  way,  the  reduction  in  propulsive 
engine  weight  in  relation  to  the  overall  vehicle  weight 
would  be  realized  without  a  corresponding  loss  in 
payload  mass  fraction  due  to  the  enhanced  effects  of  air 
drag  on  reduced-scale  vehicles. 

There  are  a  few  other  caveats  associated  with 
using  the  MEMS  combustion  results  described  in  the 
theoretical  section.  First,  MIT  demonstrated  their 
bipropellant  microthruster  with  nonhypergolic 
propellants.  For  a  TVC  microthruster  application 
where  a  small  minimum  impulse  bit  is  desired, 
hypergolic  or  catalyzed  monopropellants  would  be  the 
preferred  propellant  to  facilitate  the  intermittent  duty 
ignition  requirement  for  TVC.  Hypergolic  ignition 
followed  by  steady-state  combustion  has  not  been 
demonstrated  at  MEMS  microthruster  scales  as  of  this 
time.  The  hypergolic  or  spontaneous  ignition  delay 
would  have  to  be  accounted  for  in  the  combustion 
residence  time  analysis  discussed  previously  in  follow¬ 
up  analysis25.  This  would  involve  spark  or 
spontaneous  hypergolic  ignition  being  investigated  at 
elevated  chamber  pressures  in  order  to  achieve  steady- 
state  combustion  at  smaller  and  smaller  engine  scales. 
The  important  thing  is  that  MIT  has  demonstrated 
rocket  bipropellant  combustion  at  reduced  engine 
scales.  The  solution  to  the  spontaneous  or  spark 
ignition  problem  at  reduced  scales  is  fundamentally 
similar  in  nature  to  the  reduced-scale  combustion 
problem.  Another  caveat  is  that  MIT  demonstrated 
reduced-scale  combustion  with  gaseous  bipropellant 
injection  mixing  prior  to  combustion.  Gaseous 
bipropellant  injection  may  be  also  required  for 
hypergolic  ignition  to  be  successful.  This  brings  up 
another  important  area  for  investigation  at  reduced 
MEMS  engine  scales,  the  mixing  time  scales 
associated  with  bipropellant  liquid  to  liquid  injection 
prior  to  combustion.  Droplet  evaporation  times  should 
be  investigated  at  reduced  engine  scales  as  follow-up 
combustion  analysis. 

Thrust  to  weight  ratio  may  be  enhanced 
dramatically  at  reduced  bipropellant  engine 
characteristic  lengths,  but  the  achievement  of  a  small 
minimum  impulse  bit  as  a  result  of  reduced  engine 
characteristic  length  is  not  to  be  ignored  in  reducing 


TVC  fuel  consumption  rate  as  shown  in  Equation  29. 
The  TVC  propellant  consumption  rate,  Wp ,  was 
dramatically  reduced  for  smaller  engine  characteristic 
lengths  during  RV  RCS  attitude  keeping  dead-band 
limit-cycle  analysis.  Another  area  for  future 
investigation  would  be  regenerative  cooling  for 
intermittent  duty  microthrusters  in  order  to  maintain 

high  I SP  operation. 

As  a  final  note,  the  challenges  ahead  at 
achieving  the  bipropellant  TVC  microthrusters 
described  in  this  paper  are  formidable.  As  Prof.  Robert 
C.  Haddon  said  during  a  public  gathering  introducing 
the  new  UC  Riverside  Center  for  Nanoscale  Science 
and  Engineering26,  “To  make  contributions  to  the  field 
of  nanotechnology,  you  have  to  be  an  inch  deep  and  a 
mile  wide”.  MEMS  and  nanotechnology  involve  so 
many  disciplines,  that  the  engineer  has  to  have  a  very 
open  mind  and  to  think  outside  of  the  box  to  make 
contributions  to  microthruster  propulsion  applications. 
MEMS  and  nanotechnology  engineering  represent 
new,  and  challenging  paradigms  on  how  to  address 
traditional  propulsion  engineering27. 
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